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I. INTRODUCTION 
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We use Monte Carlo simulations to study attractive and excluded volume interactions between 
nucleosome core particles in 30 nm-chromatin fibers. The nucleosomes are treated as disk-like 
objects having an excluded volume and short range attraction modelled by a variant of the Gay- 
Berne potential. The nucleosomes are connected via bendable and twistable linker DNA in the 
crossed linker fashion. We investigate the influence of the nucleosomal excluded volume on the 
stiffness of the fiber. For parameter values that correspond to chicken erythrocyte chromatin we 
find that the persistence length is governed to a large extent by that excluded volume whereas 
the soft linker backbone elasticity plays only a minor role. We further find that internucleosomal 
^ ' attraction can induce the formation of hairpin configurations. Tension-induced opening of such 

^ I configurations into straight fibers manifests itself in a quasi-plateau in the force-extension curve 

that resembles results from recent micromanipulation experiments. Such hairpins may play a role 
\^ ' in the formation of higher order structures in chromosomes like chromonema fibers. 

PACS numbers: 87.14.Gg,87.15.Aa,87.15.La,61.41.+e 
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^H' DNA of all eucaryotic organisms is wrapped around millions of cylindrical protein spools, so-called histone octamers. 

I— I , Each complex has a radius of 5 nm and a height of 6 nm and - together with the stretch of linker DNA connecting 
to the next such spool - is called nucleosome, the basic unit of the chromatin complex [J. In a next step of DNA 
compaction the string of nucleosomes organizes itself into the chromatin fiber. For low salt concentrations a 'beads- 
on-a-string' structure is observed, sometimes referred to as the 10-nm fiber P|. For higher salt concentrations (> 40 
r^ ' niM) the fiber appears to thicken into a condensed structure with a diameter of roughly 30 nm i\ . The degree of 
f^ > compaction also depends strongly on the presence of linker histones. They cause the in- and outcoming DNA to form 
^^ a short stem- like structure '^ . In the absence of linker histones the entry-exit angle of the in- and outcoming DNA 

^P is larger, leading to more open structures. 

^^ While the structure of the nucleosome core particle (the protein spool with the 2 turns of wrapped DNA) is known 

^^ up to atomistic resolution 1^ there is still considerable controversy about the details of the structure of the 30-nm 

O chromatin fiber [BlaQilEQ- There are essentially two competing classes of models: (i) the solenoid models |3.[lOl[lll| 

' f^ and (ii) the zig-zag- or crossed-linker models jj, \12L \13^ . In the solenoid models one assumes that the successive 

I . nucleosomes form a helix with the normal vectors of the protein spools (the axis of the superhelical wrapping path 

^ ' of the DNA) being perpendicular to the solenoidal axis. The entry-exit sides of the nucleosomal DNA face inward 

^ , towards the solenoidal axis and the linker DNA must bend in order to connect neighboring nucleosomes. In contrast 

the linker DNA in crossed-linker models remains straight and connects nucleosomes that sit on opposite sides of the 

fiber. 

Chromatin fibers have been studied for various salt concentrations using electron cryo-microscopy jJ, Il4| , atomic 
force microscopy 15, .1^], neutron scattering and scanning transmission electron microscopy 17j. Structural param- 
eters such as the mass density (number of nucleosomes per 11 nm) and the linker entry-exit angle are measured to 
characterize the state of compaction. For low salt concentrations all these studies support the picture of an open 
zig-zag-like fiber structure. Similarly the force-extension curves of single chromatin fibers measured under these con- 
ditions |18l | (for a recent review see Ref. |19|) are in good agreement with the results from computer simulations 
[20l | and analytical approaches |13l l2ll |23 | based on a crossed-linker geometries of the fiber (for a recent review see 
Ref. HI). 
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The situation is less clear at physiological salt concentrations. The internal structure of the dense 30-nm fiber 
could not yet be resolved despite enormous experimental efforts, especially X-ray diffraction studies (cf. Ref. |2J| 
for a critical discussion). The interpretation of related studies on di- and trinucleosomes is also still controversial 
[a. nil I25I l2a|. Even less is known about how chromatin folds into chromosomes on larger scales. Theories for the 
elastic properties of the DNA linker backbone predict rather flexible fibers with persistence l eng ths on the order of 
10-20 nm (compared to 50 nm for uncomplexed DNA). In contrast, Langowski and coworkers |23,|23 have presented 
convincing evidence from experiments and simulations for rather large persistence length of the order of 300 nm. 
At the same time, the quasi-plateau at 5 pN that Cui and Bustamante observed in their stretching experiments at 
physiological salt concentrations points to rather delicate features in the fiber structure (for comparison, the B-S 
transition in DNA is observed at a critical force of 65 pN !29| ) . 

In the present paper we use computer simulations to study the consequence of internucleosome attractive and 
excluded volume interactions in dense crossed- linker fibers. Our chromatin fiber model is depicted on the Ihs of Fig.^ 
We model the protein spools as ellipsoidal disks and assume that linker histones induce the in- and outcoming DNA 
to form stem-like structures jj|. Our opening and rotation angles tt — 9 and are compatible with experimentally 
observed values for the 30 nm fibers mI^hI (sse the fiber "10" on the rhs of Fig. ^). As a consequence, the local fiber 
geometry remains basically unchanged when we vary the strength of the attractive interaction. Local condensation- 
decondensation transitions in chromatin fibers |l3l l23l | will be investigated in a forthcoming study J3l| . 

With respect to the degree of coarse-graining of the nucleosome structure our model is comparable to the one 
employed by Katritch et al. 20]. The main difference is their use of quenched, randomly distributed (/)-angles along 
the chain. Since we are interested in rather generic features of the dense 30 nm fiber with a narrow distribution of linker 
lengths (and hence twist angles) 1301 . we focus on the ideal case of a constant 0-angle for all linkers. Furthermore we 
follow Wedemann and Langowski [2^ in modelling the nucleosomes as ellipsoids (as opposed to spheres as in Ref. 2C|) 
interacting via a Gay-Berne potential. However, we have not included Debye-Hiickel interactions between different 
parts of the linker DNA [23 , since the Debye-Hiickel screening length is smaller than the diameter of the DNA double 
helix at physiological salt concentrations. As a consequence, intra and inter-linker electrostatic interactions can be 
accounted for by renormalized bending rigidities 32| and opening angles tt — 9 3S\ respectively. 

Problematic for all attempts to model chromatin are the soft parts of the nucleosome. The charges on the linker 
histones H1/H5 and the histone tails arc under biochemical control, allowing the cell to regulate the stem formation 
and the attractive interactions between nucleosomes |a, [sj . While simulations investigating this " tail-bridging" effect 
[35I l3q | between nucleosomes are on the way [33 , we prefer a generic nucleosome interaction- potential to a model 
which neglects the tails but accounts in detail for the surface charges on the nucleosome core [3^- Similarly, little is 
known about the stem structure (see the cryo micrographs in Ref. :4\) and for simplicity we have assumed that the 
short axis of the nucleosomal disk (the axis of the superhelical wrapping path) is oriented perpendicular to the plane 
defined by the in- and outgoing DNA linker. 

The paper is organized as follows. Section ITTl introduces our chromatin model and the methods used. We present 
our results for the mechanical fiber properties in section Hill Section Hvl provides a discussions of the observed hairpin 
structures and possible biological implications of our findings. Finally we give a conclusion in Section Ivl 

II. MODEL AND METHODS 
A. Definition of tlie chromatin model 

Our chromatin fiber model is depicted on the Ihs of Fig.^ We model the protein spools as ellipsoidal disks with 
a diameter of 10 nm and a height of 6 nm corresponding to the experimental values |5|, |39j. We assume that the in- 
and outcoming DNA are glued together in rigid stem-like structures that mimic the nucleosomal structures observed 
in the electron cryo micrographs in the presence of linker histones |j] . The stem portion is assumed to end 7 nm from 
the center of the corresponding protein spool. From this point on the entering and exiting linker DNA portions of 
length B connect to the next nucleosomes. 

The presence of a stem has to be taken into account when one calculates the length of the free linker DNA from 
the nucleosomal DNA repeat length, a number that varies not only from organism to organism but even from tissue 
to tissue of the same organism 1]. It is known that 147 base pairs (bp) are closely associated with the protein spool 
wrapping the histone core in 1.65 turns Q- We assume here that the presence of the linker histone forces the DNA 
to wrap two full turns corresponding to 177 bp. Together with the 2 nm stem this makes 189 bp that are associated 
with the core and linker histones. For instance, in the case of chicken erythrocyte chromatin the nucleosomal repeat 
length amounts to about 210 bp l] so that there are roughly 21 bp, i.e. B — 7.14 nm, of free linker DNA. In the 
following simulation runs we will always use this value for B. 

There are two angles that determine the fiber geometry: the deflection angle 6 and the rotational angle (p. The 
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FIG. 1: Lhs: Portion of the model fiber including three nucleosomes (represented in this figure as green cylinders instead of 
ellipsoids to facilitate the identification of the nucleosome orientation) connected via stems (red) to the DNA linkers (blue). 
Also indicated are the two underlying angles; the defiection angle 9 and the rotational angle (j>. Rhs: Examples of two-angle 
fibers with the arrows denoting their position in the (9, (;/>)-plane. 



former angle characterizes the entry/exit angle tt — 9 oi the linker DNA at the stem, the latter angle describes the 
rotational setting of the nucleosomal disks with respect to the DNA. Some example configurations are depicted on the 
rhs of Fig. n] with the arrows pointing to their location in the (0, 0)-space. Note that some configurations are forbidden 
since they would lead to overlapping nucleosomes; a systematic investigation of the boundary between allowed and 
forbidden structures will be provided in a forthcoming study [211 . 

In the present study we use a canonical value oi 9 = 145° which has been estimated experimentally 4] at 80 mM 
salt concentration. Having fixed B and 9 we have chosen values of that lead to a reasonable nucleosome line density. 
Based on Eq. (78) in Ref. |23 we expect 6.5 nucleosomes per 11 nm for (j) = 100° and 6.1 nucleosomes per 11 nm for 
(j) = 110° compared to experimental estimates of 6 - 7 nucleosomes per 11 nm P.[l7l|. 

The only non-rigid elements in our model are the linker DNA portions. We discretize the linkers into four segments 
in order to allow for the bending and torsional deformations. To each segment i we attach a local set of basis vectors 
{ti,ni,bi} where t^ denotes the tangent vector, n^ the normal, and b^ the binomial vector. The elastic energy of the 
linkers is then described by 
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where Ip and Itw are the DNA bending and twist persistence lengths respectively, N is the number of nucleosomes in 
our fiber and h denotes the segment length (i.e. B = 46). (3i — arccos(ti • t^+i) refers to the bending angle between 
two neighboring segments and t^ denotes the twist angle. The spontaneous bending angle psp is equal to 6 for those 
segment pairs which are connected to a stem and is zero otherwise. Tsp = (j^/^i everywhere (except at the kink where 
Tsp = 0) enforces the right-handed helicity of the DNA which in turn gives rise for the fiber twist angle (j). 

The position and orientation of the nucleosomes is calculated from the linker positions with the help of a set of 
three orthonormal basis vectors {T^, N^, Bi}. N^ is the normal vector perpendicular to the disk plane that we assume 
to be given by 



N,; = 
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Here the set {ri} denotes the vectors that connect the stems of neighboring nucleosomes (see Fig. ^. B; points from 



the tip of the stem towards the center of the nucleosomal disk, i.e. 

In -n+i] 

FinaUy, T,; is defined as 

T, = N, X B,. (4) 

For the interactions of the nucleosome core particles we use the same variant of the Gay-Berne (GB) potential for 
ellipsoids of arbitrary shape !4Qi, 41j| as in a recent study of a stacked-ellipsoid model of DNA _42j . Here we set the 
structure matrix (Eq. (15) in Ref. ji^l) equal to 

^h/2 \ 
5 = ( d/2 , (5) 

d/2/ 

where d = 1.67h is the diameter of the nucleosome core particle with a canonical value for the core particle height of 
h — 6 n m p . The effective distance of closest approach hi2 between two ellipsoids is calculated using Eqs. (17)-(19) 
in Ref. J43. In cases where we study purely repulsive hard core interactions, we reject all Monte Carlo moves leading 
to values /112 < 0. In cases where we study attractive interactions, the distance dependent part of the interaction 
potential Eq. (14) in Ref. [ij is given by a modified Lennard- Jones potential 

The orientation dependent parts are defined in Eqs. (20)- (25) in Ref. [431 . This parameterization of the GB potential 
leads to a preferred lateral spacing of about 7.0 nm and a vertical spacing of about 11.0 nm of the disks. The values 
are in good agreement with experimental results for columnar phases of nucleosome core particles |22| • In this paper 
we characterize the strength of the nucleosome-nucleosome attraction by the depth of the attractive well for two 
parallel disks e = lAlecB instead of the prefactor ecB m Eq. (0. 

B. Methods 

We use a Monte-Carlo scheme to simulate the chromatin fiber which relies on the following three moves: (i) a local 
move where one randomly chooses a nucleosome that is rotated around an axis determined by two points on the in- 
and out-coming linker DNA by a small random angle, (ii) a non-local pivot move where a random segment point is 
chosen at which the shorter part of the chain is rotated around a random axis by a random angle and (iii) a non-local 
crankshaft move where two random points along the DNA segments define the axis of rotation around which the inner 
part of the chain is rotated. The moves are accepted or rejected according to the Metropolis scheme p^ . 

We simulated either fibers with A^ = 50 (simulations with applied stretching force) or A^ = 100 nucleosomes 
corresponding to 10 and 20 kbp of DNA respectively. The number of nucleosomes in pulling experiments by Cui 
and Bustamante ^^ is on the same order (w 300 nucleosomes). Each simulation run consists of 200000 MC sweeps, 
where one sweep corresponds to Njjna trials with Njjna being the number of DNA segments. The amplitudes are 
chosen such that the acceptance rate equals approximately 50%. Every 20 sweeps we save a configuration. As initial 
conformation we used the relaxed (T = 0) fiber structure. In order to determine the longest relaxation time Tcorr 
of the system, we measured the 'time' correlation functions of the energy, the mass density, the end-to-end distance 
and the twist. We typically find that Tcorr ~ 100 MC sweeps. Note, however, that strong attractive interactions can 
cause a glasslike trapping of the fiber in local energy minima. As a consequence, only a limited range of attractive 
well depths can be studied meaningfully in simulations. In the present case, the simulation runs for e = 4 ksT were 
clearly non-ergodic for low stretching forces. In this case we averaged the results over 10 independent simulation runs. 

The calculation of the elastic constants of the simulated fiber is performed as follows. First we determine the fiber 
axis by calculating the centers of mass Ci = ll/Nc) Yl^j=i " ^j ^^ groups of Nc neighboring nucleosomes. Nc is chosen 
to match approximately one or two helical turns. Then we calculate the autocorrelation function (t^ -tj) of the tangent 
vectors of the fiber axis 

t. = ^±^ (7) 
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FIG. 2: Effect of the excluded volume interaction on the bending persistence length of the fiber for two deflection angles 9 with 
(j> = 100° and B — 7.14: nm. For a given linker length B the fiber persistence length Ip [nm] grows with increasing disk size. 
For very small disk sizes, the measured persistence lengths converge to the analytical result, Eq. (97) of Ref. [2^. On the rhs 
we show some snapshots of simulated fibers with varying disk sizes for 6 = 145° (top row) and — 90° (bottom row). 



and extract the persistence length of the fiber Ip from an exponential fit 
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The stretching modulus kgT^ can be estimated via 
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with AL = L — (L) being the mean deviation from the average contour length of the fiber (defined as the length of 
the fiber axis L = X]i=i " l^i+i — Ci\). It should be noted that depending on Nc the estimated values of the contour 
length, the persistence length and the stretching modulus vary. In case of Nc being too large the stretching modulus 
is underestimated since bending fluctuations within i and i + Nc are averaged out so that the contour length of the 
fiber appears to be smaller. On the other hand, values of Nc which are too small lead to a helicoidal fiber axis, and 
the contour length of the fiber is overestimated. This entails a systematic error which must be minimized. We found 
that values of Nc corresponding to one or two helical turns lead to reasonable estimates. 



III. RESULTS 



In order to quantify the influence of the nucleosome excluded volume on the flber stiffness we performed a series 
of simulations with purely repulsive interactions between the nucleosome core particles, constant linker geometry and 
variable core particle size. The Ihs of Fig. [3 shows our results for the persistence length for two fiber geometries; the 
corresponding snapshots in the rhs of the figure illustrate the observed fiber stiffening with increasing core particle 
size. 

To be more specific, we studied two different sets of angles, namely 9 — 145° and (p — 100° as well as 9 = 90° 
and (j) = 100°. In both cases, starting at ft, = (i.e., no nucleosomes present) we observe the theoretically expected 
values for the persistence length of the linker backbone (Eq. (97) in Ref. ^i), namely Ip ~ 14 nm and Ip ~ 34 
nm, respectively. We verified this formula also for many other values of 9 and (j) (data not shown). However, with 
increasing h the value of Ip increases and reaches values of 240 nm and 90 nm, respectively, at the canonical value of 
h = 6 nm. This corresponds to a 17-fold and 3-fold increase respectively of the fiber stiffness relative to the theoretical 
prediction. The different degrees of stiffness (and stiffening) directly reflect the different nucleosome densities in the 
two cases: in the 9 = 145°-fiber the nucleosomes are always quite close to each other whereas the 9 = 90°-fiber is still 
much more open. 
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FIG. 3: Force-extension curves for fibers with 9 = 145°, </!) = 110° and B = 7.14 nm. The curves correspond to fibers with 
different nucleosomal attraction, namely e = (pure hardcore) and e = 2, 3, 'iksT. At e = AksT occurs a force plateau around 
2pN . Snapshots of fibers with 6 = 145°, = 110°, and B — 7.14 nm at different stretching forces for e = AksT. To facilitate 
the detection of hairpins one half of the chain is shown in green while the other half is shown in cyan. The end nucleosomes 
are labelled red. The fiber at / = pN shows a kink close to the center of the chain. Up to / = 2 pN the kink is still present 
but the ends get pulled out. For / = 15 pN the fiber is stretched and nucleosomal contacts are broken. 

We also performed simulations of stretched fibers with 9 — 145° and </> = 110°. The resulting force-elongation 
curves (red symbols in Fig. ^ and histograms of the end-to-end distance distribution (red symbols in Fig. ^ show 
the typical behavior of extensible worm-like chains. The entropic small force regime is strongly suppressed, since 
the / = contour length, L = 85nm, of fibers containing 50 nucleosomes is much smaller than their persistence 
length of Ip — 220 nm. The effective stretching modulus 7 = 8nm~^ at finite extensions is smaller than the value 
7 = 14nm^^ deduced from the analysis of the length fiuctuations (see Eq. O and larger than the the theoretical 
prediction 7 ~ 3nm^^ which is based on the linker mechanics alone. 

For our study of the effect of attractive interactions between core particles we focused on linker backbone geometries 
which by themselves already lead to relatively dense fibers {9 = 145° and (p = 110°). As a consequence, the local fiber 
geometry remains basically unchanged when we vary the strength of the attractive interaction from e = (the case 
of purely repulsive interactions discussed above) to e = 4. For example, the observed attraction induced reductions 
in the fiber contour lengths are only of the order of 10%. 

Weak attraction up to e = iksT (green symbols in Figs. |31 and ^ has only a small effect on the observed fiber 
properties. While the contour length decreases to about L = 80nm, there is a corresponding small increase in both 
the fiber stiffness and its stretching modulus, yielding an overall small increase in the end-to-end distance. The 
situation changes dramatically for larger values of the attractive well depth. In the absence of external forces fibers 
with e = iksT and, in particular, with e = AksT (blue/magenta symbols in Figs.Eland^ have considerably smaller 
and more broadly distributed end-to-end distances Re- For small forces the e = '^ksT fiber shows a quasi-plateau in 
the force-extension curve. 

The origin of this peculiar behavior can be identified by inspecting snapshots of the e = 'iksT fiber (Fig. O: At 
f — pN the fiber forms a hairpin structure where the end nucleosomes (in red) are located at the same end of the 
structure. The hairpin persists at / = 0.5 pN but opens up into a straight configuration around / = 2 pN. For larger 
forces all observed force elongation curves show qualitatively the same WLC behavior. 



IV. DISCUSSION 



Our results show that the properties of dense chromatin fibers are dominated by the interactions between the 
nucleosome core particles, while the mechanical properties of the linker backbone play a less important role. The 
dominant effect of excluded volume interactions is to stiffen the fiber. The persistence lengths of the order of 200 - 
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FIG. 4: Probability density of the end-to-end distance Re of a fiber with = 145°, = 110°, and B 
energy well depths e and stretching forces / as specified in the legends. See text for details. 



7.14 nm for various GB 



250 nm observed in the present study are in good agreement with the values obtained by Langowski et al. [23, |28| . 
The surprising result of our simulations is that attractive interactions of a few fesT per nucleosome pair are sufficient 
to bend fibers with a contour length of about 1/3 of the persistence length into dense hairpin configurations. Judging 
from Fig. 5 in Ref. ;20] hairpins also occurred in the simulations by Katritch et al. - even though the authors did 
not discuss this issue in their paper. Since the observed force-elongation curves show qualitatively similar features 
in all three cases^it is tempting to speculate that similar effects also occurred in the stretching experiment by Cui 
and Bustamante jl^l- In the following, we will discuss the condensation of semi- flexible filaments, the local structure 
of the hairpins observed in our simulations, the signature of hairpin opening in stretching experiments, and possible 
implications of our observations for the folding of chromatin into chromosomes. 



Condensation of semi-flexible filaments 



The condensation of semi-flexible filaments due to short-range attractive interactions has recently been treated in 
detail by Schnurr et al. |44| . They introduced a condensation length Lc = ^/ksTlp/aattr (o'attr'- attraction energy 
per length) by balancing expressions for typical bending and surface energies, ksTlp/L ~ UattrL. Filaments with a 
contour length L/ Lc < 0(10) remain extended, while longer chains aggregate into structures with a typical size of 
the order of Lc- So-called "racquet" states consist of a straight stem where the chain is folded back on itself several 
times and where 180° U-turns are shaped like the head of a tennis racquet. A racquet state with one turning point 
resembles our hairpin structure. Longer fibers lower their energy by having multiple turning points. For infinitely 
thin filaments Schnurr et al. were able to show that the kinetically preferred racquet states |43 have slightly higher 
energies than toroidal structures. 

In the present case, we measured a typical attraction energy per length of aattr ~ {3/4) kBTnm^^ for hairpin 
conformations of the e — AkBT-Ghei in good agreement with an estimate of e/h = AksT/Gnm for the attractive 
energy per nucleosome pair of height h. Using this estimate and a persistence length Ip = 240nm, chromatin fibers 
with weak attractive interactions should start to condense, if they are longer than about C'(10)ic ~ ^^fesT/e 400 nm. 
Note that the above argument breaks down for e > GfcsT, when the predicted racquet head or torus radii become 
smaller than the radius of the 30 nm- fiber. 




FIG. 5: Contact matrix of the fiber with e = 4 ksT, 9 = 145°, (j> = f 10°, and B = 7.14 nm. (a) At an external tension / = 3.5 
pN the fiber is stretched and there is no hairpin, (b) to (d) Without an external force, / = pN, one can clearly detect 
hairpin structures. The simulation runs (b) and (c) show the occurrence of two hairpins whereas in (d) only one hairpin can 
be identified. 

Can the structure of the hairpins we observe really be understood in the framework of the WLC model? Or does 
one have to consider the relatively sharp turning points with the curvature radii of the order of the fiber radius as 
localized kink defects? Following the analysis of Schnurr et al. |4J|, the fibers in our simulations are too short to 
fold back on themselves. The WLC estimate of the bending energy for a U-turn with radius 15nm is Ukink ~ 25kBT, 
while the onset of hairpin formation at e = iksT in 80 nm long fibers suggest that the actual energy penalty for 
the kinks is as small as Ukink ~ Ofc^T. We note that the combination of finite fiber radii and relatively small kink 
energies should stabilize racquet structures relative to toroids. 



B. Hairpin structure 



The internal structure of hairpin configurations shows up very clearly in the contact matrix A4 contact of the fiber that 
we depict in Fig.[Sl For a given configuration the contact matrix is defined as follows: If a pair of nucleosomes i and j 
is in contact (center-to-center distance smaller than 1.5d = 15nm) one has Mcontact = 1; otherwise M-cgntcict = 0- By 
taking the average of the matrix elements for many configurations of a given simulation run we obtain two-dimensional 
histograms as it is shown in Fig. |S| For large forces no kinks are present and the nucleosomes of the fiber form a 
quite regular crossed-linker structure. We find two pronounced diagonals at the positions i — j ±2 that correspond 
to "short-ranged" excluded volume interaction. Less pronounced side stripes parallel to the main diagonal are also 
observed and describe interactions between nucleosomes of neighboring turns, especially at i = j ± 5 and at i = j ± 7. 
Hairpins manifest themselves by crosslike patterns in the contact matrices with branches along the secondary diagonal. 
The location of the kink corresponds to the point where the two diagonals cross each other. Closer inspection of the 
histograms shows that the hairpin arms outside the kink region are partially interdigitated but not strongly perturbed. 



However the effect is strong enough to effectively prevent any possibility of sliding of the two arms with respect to 
each other on the time scale of our simulations. Therefore a hairpin structure, once formed, can be very stable with a 
frozen in position of the turning point. With respect to the fiber structure at the turning point, Fig.E^d) in particular 
shows disruptions of the internal order of the nucleosomes: The diagonals i — j ± 2 are broken at that position, 
presumably due to broken contacts of nucleosomes located at the outside of the turning point. This observation 
supports the picture of a localized kink defect. 

C. Large scale vs. local decondensation in stretching experiments 

Toroidal or racquet-like condensates become unstable under the influence of a stretching force, if the applied force 
/ exceeds a value on the order of the attractive energy per unit length, aattr ■ Neglecting kinetic barriers (for a 
discussion of DNA unwrapping from nucleosomes or toroidal aggregates see H^), the corresponding force plateau in 
the force-extension curves for chromatin fibers should be observed for forces of / « (Jattr ~ (2/3)(e/fcsT)pN. In order 
to experimentally distinguish this scenario from local structural changes which are expected to occur at similar force 
levels |l3| . one should plot force- relative extension curves for fibers reconstituted under similar conditions on DNA 
strands of different length. In the case of local structural changes, the obtained curves coincide, while for a global 
decondensation the size of the unperturbed state is of order Lc independently of the fiber length. 

D. Biological implications 

As a final point we discuss possible biological implications of fiber stiffening and hairpin formation. Both effects are 
associated with characteristic length scales; fiber stiffening with a persistence length on the order of 300 nm, hairpin 
formation with a condensation length on the order of a few tens of nm. 

We note again that the biochemical control of the effective interaction between nucleosomes is essential for the 
condensation and decondensation of chromosomes during the cell cycle |43| • There is indirect experimental evidence 
for the formation of chromosomal DNA loops having a size of about 50 kilobases [ig (corresponding to about 250 
nucleosomes), e.g. from the comparison of separation patterns of excised large DNA fragments by pulsed field 
elec trop horesis with the patterns obtained by DNA cleavage by topoisomerase II of the nuclear matrix of the cell 
p9ll50| . Using a nucleosome density of 6/(11 nm), 250 nucleosomes correspond to a fiber length of 460 nm, i.e. about 
twice the persistence length of 240 nm found in our simulations. 

There are two scenarios linking the persistence length to the loop size; both are based on the observation that 
the free energy cost for the formation of a loop is smallest for semi-flexible chains with a length of about twice the 
persistence length. For shorter chains the bending energy is high, for larger chains the entropy of chain conformations 
counteracts loops formation. In the first (equilibrium) scenario, the chromatin fiber loops around strongly attracting 
organizing centers. Spherical centers could induce the formation of rosette structures [5l|. Filamentous organizing 
centers could lead to structures resembling cartoons of chromosomes often found in biological textbooks, cf. e.g. 
Ref. |53| . The second scenario which links the persistence length to the loop size is closer to what we observed in our 
simulations: if the attraction between fibers is strong enough to prevent structural reorganization after the two halves 
of the fiber have tightly closed the gap between the point of first contact and the hairpin defect, then loops with the 
typical size of the fiber persistence length are preferred for kinetic reasons. Hairpins have indeed been observed in 
cryo-EM pictures of chromatin fibers in the presence of MENT [531, a protein that is involved in the formation of 
dense, transcriptionally inert sections of chromatin, the so-called heterochromatin. 

If, on the other hand, structural reorganization is still possible in collapsed fibers, then the size of the aggregates 
should be controlled by the smaller condensation length and not by the persistence length of the 30-nm fiber. Fibers 
with diameters ranging between 60 to 130 nm diameter have indeed been observed in mitotic and in Gl chromosomes 
and are called chromonema fibers |5J|, cf. also Ref. [Ig. However, the folding of 30-nm fibers into superfibers with 
a diameter on the order of the condensation length requires a mechanism that prevents the formation of globular 
aggregates. 

V. CONCLUSION 

We have used computer simulations of model chromatin fibers to investigate the influence of excluded volume and 
attractive interactions between the nucleosomes on the large-scale structure and elasticity of the 30 nm chromatin 
fiber. Our results shed light (i) on the discrepancy between the theoretically expected and the observed persistence 
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length of 200-300 nm as well as (zz) on the somewhat counterintuitive observation that such fibers nevertheless seem 
to be able to curl up into chromonema fibers with diameters of 60 - 130 nm. 

Our results clearly show that the stiffness of dense fibers is dominated by excluded volume interactions between 
nuclcosomes. The observed persistence lengths exceed estimates based on the linker backbone elasticity by one order 
of magnitude. With respect to internucleosomal attraction we have concentrated on generic, non-local aspects of 
fiber condensation. The 2 pN force plateau observed in our simulations reflects a structural feature beyond the 30 
nm fiber: the opening of a hairpin. Our results suggest that aligned (anti-) parallel fibers arc only weakly perturbed, 
but that the fiber geometry is locally disrupted at the 180° turning point. The crucial point is that the energetic 
cost of this localized defect is considerably smaller than for a smoothly bending semi-flexible filament with the same 
bending persistence length as chromatin. This effect could be important for an understanding of chromatin folding 
in chromosomes. 
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